Mon. Not. R. Astron. Soc. 000, 



ism 



Printed February 1, 2008 (MN WT^K style file v2.2) 



Chandra observations of the galaxy cluster Abell 1835 



R. W. Schmidt,* S. W. Allen and A. C. Fabian 

Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CBS OH A, United Kingdom 



o 
o 



Received 



o 



43 
6 

u 

CO 

C3 



ABSTRACT 

We present the analysis of 30 ksec of Chandra observations of the galaxy cluster Abell 
1835. Overall, the X-ray image shows a relaxed morphology, although we detect sub- 
structure in in the inner 30kpc radius. Spectral analysis shows a steep drop in the 
X-ray gas temperature from ~12keV in the outer regions of the cluster to ~4keV in 
the core. The Chandra data provide tight con straints on the gravitational pot ential of 
the cluster which can be parameterized by a Navarro, Frenk & White (1997) model. 
The X-ray data allow us to measure the X-ray gas mass fraction as a function of radius, 
leading to a determination of the cosmic matter density of f2 m = 0.40±0.09 fcjo ■ The 
projected mass within a radius of ^150 kpc implied by the presence of gravitationally 
lensed arcs in the cluster is in good agreement with the mass models preferred by 
the Chandra data. We find a radiative cooling time of the X-ray gas in the centre of 
Abell 1835 of about 3xl0 8 yr. Cooling flow model fits to the Chandra spectrum and 
a dcprojcction analysis of the Chandra image both indicate the presence of a young 
cooling flow (~ 6xl0 8 yr) with an integrated mass deposition rate of 230^5qM q yr^ 1 
within a radius of 30kpc. We discuss the implications of our results in the light of 
recent RGS observations of Abell 1835 with XMM-Newton. 

Key words: galaxies: clusters:general - galaxies: clusters: individualAbell 1835 - 
cooling flows - intergalactic medium - gravitational lensing - X-rays: galaxies 
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1 INTRODUCTION 

Since it s launch in 1999 July , the Chandra X-ray obser- 
vatory |(Weisskopf el til. 200C) has ubseived a laige iium- 
ber of 



galaxy — clusters, many — of which are amongst the 
most massive objects of their kind. In particular, high res- 
olution 



spalially-iesulved speeliuseupy wiLli the Advanced 



CCD I maging Spectrometer (ACIS) on Chandra is leading 



to significant advances in the way we uiideislaiid these sys- 
tems. 



As 



lliislialed in a recent papei by Allen el al. (2001a ) 



Chan dr a observations enable U3 to constrain the proper — 
ties of galaxy clusters with unprecedented accuracy. Chan- 
dra data permit direct measurements of the X-ray gas den- 
sity, temperature and total mass profiles at a resolution of 
~ 20 hj kpc for the most massive clusters observed at red- 
shifts z ~ 0.2, paving the way for detailed cosmological stud- 
ies. 



est Cluster Sample ( Ebeling et al. 1998 ) and has previously 
been inferred to contain a strong cooling flow (with a nom- 
inal mass deposition rate of ~ 1000 Mq yr~ ; e.g., Allen 



et al. 1996). The central dominant galaxy exhibits powerful 



optical emission lines and a strong UV/blue continuum asso- 
ciated with large amounts of ongoing star formation (Allen 



1995 



Crawford et al. 1999). The observed Balmer em ission 



line r atio also indicates significant intrinsic reddening ( Allen 
1995|, E(B-V) = 0.49JI°;^). This is consistent with the de- 



tection of 850 /im emission from the cD galaxy by Edge et 



al. (1999), which those authors attribute to emission from 



In 



warm dust heated by the young, hot stars. 

The great concentration of mass in the cluster is also 
revealed by several features in optical images that can be 
identified as gravitationally lensed background objects. In 
particular, a large arc was discovered by Edge e t al. (in 
preparati on) to the south-east of the cD galaxy. Allen et 



this paper ™. present Qumdm observations of the aL ( 1996 D used this . arc to P ut constraints on the cluster 



galaxy cluster Abell 1835 at a redshift z = 0.2523. This clus- 
ter is the most luminous system (Lx = 2.7 x 10 45 ergss _1 
in the 0.1-2.4 keV ROSAT band) in the ROSAT Bright- 
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mass inside the radius of the arc. We will discuss several 
further possibly lensed objects in this study. 

The structure of this paper is as follows: Sect. de- 
scribes the Chandra observations and Sect. ^ the basic imag- 
ing results. Sect. W discusses the spatially-resolved spec- 
troscopy of the cluster using the Chandra data. Sect. H 
presents a detailed mass analysis using the independent X- 
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ray and lensing data.In Sect. H we examine the properties 
of the central cooling flow and compare our results with 
those from recent observations made with the XMM-Newton 
X-ray observatory. We present our conclusions in Sect, bj 
All quantities are given throughout using a cosmology with 
Ho = 50kms _1 Mpc -1 and qo = 5. Unless otherwise noted 
all error bars are 1 a (68.3%) confidence intervals. 



2 OBSERVATIONS 

Chandra observed Abell 1835 on 12 December 1999 for a 
total of 19.6 ksec. The target was centred in the middle of 
node 1 of the back-illuminated ACIS chip S3 (ACIS-S3), 
near the nominal aim point of this detector. The detector 

o 

temperature during the observation was -110 C. A further 
10.7 ksec Chandra observation of Abell 1835 was obtained 

o 

on 29 April 2000 at a detector temperature of -120 C. In 
this second exposure, the cluster was centred close to the 
node boundary in the geometrical centre of the ACIS-S3 
chip, in the middle of a group of five dead columns. This 
has prevented us from using this second data set for the 
spectral analysis presented in this paper. We have, however, 
used both data sets to produce the images in Figs, y and ^ 



3 IMAGING 

In Fig. Ill a contour plot of the 30 ksec Chandra exposure 
of the inner region of the galaxy cluster in the 0.3-7.0 keV 
energy band is shown with a side length of 3.3 arcmin- 
utes (IMpc at the distance of the galaxy cluster). About 
80,000 photons are collected in this image. For comparison, 
an archival 20min B-band exposure taken on 27 February 
1998 with the Canada France Hawaii Telescope (CFHT) of 
the central 1.5' x 1.5' is shown in Fig. H. In the optical image, 
the large arc (object A) discovered by Edge et al. (in prepa- 
ration) can be seen to the south-east of the central cluster 
galaxy. This arc is thought to be a gravitationally lensed 
image of a galaxy far behind Abell 1835. The presence of 
this arc, notwithstanding the large X-ray emissi on, suggests 



that A bell 1835 is a very massive galaxy cluster. Allen et al 



(1996) estimated the projected mass inside the radius of the 
arc to be between 1.4xl0 14 M and 2xl0 14 M Q , if the arc 
has a redshift greater than 0.7 and is situated on a critical 
line of the gravitational potential of the cluster. In Fig. H 
a smaller arc can be seen to the south-west (object B) of 
the central cluster galaxy. The arrow to the north-west of 
the cD galaxy points to a radial feature which we discuss in 
Sect. |J. 

1.5 arcmin to the 



Except for an apparent point source < 



cast (this i s discussed in Fabian et al. 2000a and Crawford 
ct al. 200C), the Chandra image is dominated by emission 



from the hot intracluster gas of the galaxy cluster. Overall, 
the cluster has a regular, relaxed morphology. The isophotes 
outside the core are elliptical with an axis ratio of minor 
and major axis ~ 0.85. In Fig. CH the average counts per sec 
per arcsec 2 as a function of radius are shown as counted in 
annuli around the galaxy cluster centre in the 0.3-7.0 keV 
band. The profile shows the steep surf ace brightnes s profile 
typical of cooling flow clusters (e.g., Fabian 1994). Note, 




Figure 1. X-ray contours of the central 3.3x3.3 arcmin (lxl Mpc) 
of Abell 1835 in the 0.3-7.0 keV band. The contours were drawn 
at 1, 2, 4, 8, 16, 32, 64, 91 and 108 counts per lxl arcsec pixel. 
The coordinates are J2000. 




Figure 2. Optical B-band image of the inner 1.5' X 1.5' 
(0.45x0.45 Mpc) of Abell 1835. North is up and east is to the 
left. The cross marks the galaxy core. The greyscale has been 
wrapped once for clarity. 



however, that Chandra detects a flattening of the profile 
within r ~ 30kpc. 

In Fig. ^ the central 30"x30" (150 kpc side length) of 
Fig. hi are shown. The image was smoothed with a Gaussian 
kernel with a full-width at half-maximum of 1.6 pixels (0.8 
arcsec). In the inner 20 arcsec (100 kpc) of the cluster Chan- 
dra resolves substructure in the X-ray emission. There are 
a couple of main features visible: a well-defined peak just 
east of the geometrical centre (denoted by the black ring; 
J2000.0 coordinates given by Chandra: RA 14:01:02.07, DEC 
+2:52:43.2) of the cluster, a more inhomogeneous structure 
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^ radius (kpc) 

Figure 3. Surface brightness profile of Abell 1835 as a function 
of projected separation from the cluster centre. The count rates 
have been azimuthally averaged in annuli of width of 1 arcsec 
(4.9 kpc). The data points are positioned on the abscissa at the 
average of the inner and outer radius of the annuli. The data 
points are shown with Poisson error bars. 



to the south-west of the geometrical centre, between these 
two regions an elongated strip with a reduced count rate, 
and a small emission front to the south-east of the cluster 
core (between the arrows in Fig. g). The brightness gradient 
towards the north-west is steeper than towards the south- 
east. These features may be the sign of a merging subclump. 
The number of counts per half-arcsecond pixel in the bright- 
est regions in the combined 30ksec of our data is between 
30 and 40. This means that the detailed structures need to 
be analysed with care because the Poisson error is still of 
the order of 20%. The main features, however, are robustly 
detected and indicate the presence of variations of the X-ray 
emission or absorption in the core of Abell 1835 on scales 
of several kiloparsecs. Comparison with the optical image 
(Fig. 0) shows that the emission inhomogeneities are situ- 
ated inside the central cluster galaxy. 



4 SPATIALLY RESOLVED SPECTRAL 
ANALYSIS 

4.1 Data extraction 

We extract spectra from selected regions of the longer 19.6 
ksec exposure Chandra data set using the CIAO software 
package distributed by the Chandra X-ray Observatory Cen- 
tre (CXC, web site: cxc.harvard.edu). 

We only work with photons detected in the ACIS-S3 
chip, which is also known as ACIS chip 7. This is a back- 
illuminated chip and does not suffer from the degrada- 
tion of energy resolution that affects the front-illuminated 
ACIS chips on the Chandra observatory. All significant point 
sources were masked out and excluded from the analysis. 

Response matrix files (RMF) and ancillary response 
files (ARF) were produced using a PERL script generously 
provided to us by Roderick Johnstone. This script combines 
individual RMFs and ARFs for each 32x32 pixel region on 
the ACIS-S3 chip (generated using the mkrmf tool in the 
CIAO software suite available from the CXC website) in a 
photon number-weighted manner. 

Using the task grppha from the FTOOLS software 
package provided by the NASA High Energy Astrophysics 
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Figure 4. X-ray image of the central 30" x30" (150kpcxl50kpc) 
of Abell 1835 in the 0.3-7.0 keV band. North is up and east 
is to the left. The image was smoothed with a Gaussian with 
a full-width at half-maximum of 0.8 arcsec. The pixels are 
0.5x0.5 arcsec. Only pixels with at least 6 counts are plotted. The 
circle marks the geometrical centre of the cluster mentioned in the 
text. The cross marks the absolute posit ion of the galaxy core (see 
Fig. ph determined by Allen et al. (1999) when plotted in the coor- 
dinates given by Chandra. The uncertainties associated with the 
galaxy position and the Chandra astrometry are ~ 1 arcsec. The 
two arrows point at the small emission front extending between 
them. The greyscale has been wrapped once for clarity. 



Science Archive Research Center (HEASARC, FTOOLS 
web site: heasarc.gsfc.nasa.gov/lheasoft/ftools/) we have 
grouped the extracted spectra so that we have at least 20 
counts in a bin corresponding to a certain range of energy 
channels. We can then use Poisson error bars to account for 
the uncertainty in the number of photons per bin. 

We have have generated background spectra using 
Maxim Markevitch's method and program that is available 
from the CXC web site. We extract the background data 
from the same regions on the chip as the spectra extracted 
from the Abell 1835 data set. We also scanned the light 
curve of the Abell 1835 data for flares or background activ- 
ity using the tool lc_clean provided by Markevitch. The light 
curve is stable throughout the 19.6 ksec exposure. 



4.2 Colour profile analysis 

In order to illustrate the spectral energy distributions from 
different regions in Abell 1835 we have firstly constructed 
an X-ray colour profile. Two separate images were created 
in the energy bands 0.5 — 1.3 and 1.3 — 7.0 keV (0.6 — 1.6 
and 1.6 — 8.8 keV in the rest frame of the source) with a 
1.97 arcsec (4 raw detector pixels) pixel scale. These soft 
and hard X-ray images were background subtracted and flat 
fielded (using the exposure map tools available on the CXC 
web site and taking full account of the spectral energy dis- 
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Figure 5. X-ray colour profile of the central 1.7arcmin (500 kpc) 
of Abell 1835. In this figure the ratio of the azimuthally averaged 
counts in the 0.5 — 1.3 and 1.3 — 7.0 keV bands (0.6 — 1.6 and 
1.6 — 8.8 keV in the rest frame of the source) is plotted as a 
function of the radius of the annuli. 
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Figure 6. Temperature kT as a function of radius r from the 
cluster centre out to 3.3 arcmin (IMpc). For data points at radii 
Js 30arcsec (150 kpc), we have plotted the temperatures derived 
using models with the absorption fixed at the Galactic value. 



tributions of the detected photons). Azimuthally-averaged 
surface-brightness profiles for the cluster were then con- 
structed in each energy band, centred on the geometrical 
centre of the X-ray emission. The colour profile formed from 
the ratio of the surface brightness profiles in the soft and 
hard bands is shown in Fig. pi 

From examination of Fig. H we see that at large radii 
the observed X-ray colour ratio is approximately constant, 
with a mean value of 1.05 ± 0.03 (la error determined from 
a fit to the data between radii of 200 — 400 kpc). Within 
a radius of ~30 arcsec (150 kpc), however, the colour ratio 
rises, indicating the presence of cooler gas. 



4.3 Single-temperature model 

We investigate now the temperature profile that causes the 
colour profile shown in Fig. H. To do this we have extracted 
spectra from annuli around the geometrical centre (marked 
by the black circle in Fig. Mj) of the X-ray emission with 
photon energies in the 0.5-7.0 keV band. We use a slightly 
smaller energy band for the spectral analysis than we used 
for the images in order to get robust answers at soft energies. 
The annuli have widths of 5, 10, 25, 50 and 100 arcsec (25, 
50, 125, 250 and 500 kpc, respectively) and provide sufficient 
counts (see Tab. \w to study the spectral properties of the 
X-ray gas. 



We use the software package XSPEC ( Arnaud 1996 ) to 
model the spectra with the MEKAL plasma emission model 
([Kaastra fc Mewe 1993J; JLicdahl et al. 1995|) multipli ed by 



the PHABS (Balucinska-Church & McCammon 1992) pho- 



toelectric absorption model 

MODELi = PHABS (JV H ) X MEKAL (T; Z; K). 



The free model parameters are indicated in brackets: the 
absorbing equivalent hydrogen column density iVn, the gas 
temperature T, the metallicity Z (in units of solar abun- 
dance), and a normalization factor K. 

An account of the relevant parameters and the model 
fits is given in Table h|. All our fits are statistically accept- 
able, with reduced \ values of around or slightly above 
unity. Since the absorbing column density and the metallic- 
ity are not well constrained in the outer annuli, we have also 
carried out a second analysis in which we grouped the data 
into three larger regions. This procedure eliminates free pa- 
rameters and thus yields more robust estimates. The results 
are shown in Table ti. 

It can be seen from Table that the spectra from annuli 
with inner radii larger than 30 arcsec (150 kpc) are consis - 
tent with Galactic absorption (Dickey & Lockman 1990). 



(1) 



We have thus repeated the isothermal model fits according 
to eq. dll) for these annuli with fixed Galactic absorption. 
The corresponding best-fit models are also given in Table III 
For the outermost annuli, this leads to a smaller predicted 
temperature with smaller uncertainties. 

In Fig. m the results from the model fitting for the am- 
bient gas temperature T are plotted. For annuli with inner 
radii ^30 arcsec (150 kpc), we have fixed the absorbing col- 
umn density to the Galactic value. The best-fit temperatures 
reveal a strong drop from kT ~ 12 keV in the outer regions 
of the cluster down to 4 keV in the centre of the cluster. We 
note that the photons in the annuli have been emitted in a 
hollow cylinder with a cross-section given by the annulus. 
The measured temperatures are fits to the spectrum from 
this whole region and thus have to be regarded as average 
(emission weighted) temperatures. 

We have plotted the metallicity and absorption re- 
sults of Table H in Fig. u\ The metallicity appears approxi- 
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Table 1. Isothermal model fits. The model described in eq. (hi) with the free parameters equivalent hydrogen column density Nn, 
temperature T and metallicity Z was fitted to spectra extracted from annuli around the geometrical centre of Abell 1835 with inner 
radii r% and outer radii T2- For annuli with ri Js 250 kpc (50arcsec) we have modelled the X-ray spec tra from the annuli both with free 
equivalent hydrogen column density Nn (top block of models), and with fixed Galactic absorption (Dickey & Lockman 199C, bottom 



block of models). Since neither the metallicity nor the absorbing column density can be smaller than zero, the lower limit has always been 
fixed accordingly for these parameters. In columns 6 and 7 we list the total number of counts and the expected number of background 
counts in the 0.5-7.0 keV band. In the last three columns the x 2 values, the number of degrees of freedom (DOF) and the reduced 



X 



X 2 /DOF are given. 



r\ (kpc) 


r-2 (kpc) 


kT (keV) 


Z (solar) N 


H (10 20 cm- 2 ) 


counts 


exp. background 


x 2 


DOF 


x 2 





25 


4 o+ - 3 
^■ u -0.3 


24+0- 07 
u,z4 -0.07 


Q QO+0.82 

°- yo -0.82 


4791 


2.1 


148.1 


122 


1.21 


25 


50 


4 7+0-2 
**'-Q.a 


0Q1+0.05 
u - ol -0.06 


o Q7+0.66 
°- a '-0.64 


8300 


8.5 


218.1 


161 


1.35 


50 


100 


7 k+0-5 
'•°-0.4 


n 9 o+0.06 


o 97+O.58 

■^'-o.ss 


10776 


27.0 


233.8 


193 


1.21 


100 


150 


8 5+ - 7 
°-°-0.7 


3Q+OIO 

u -° 3 -o.io 


9 OQ+0.66 
z - z3 -0.65 


7680 


51.2 


186.3 


172 


1.08 


150 


250 


Q 1+0-8 

^■'--o.r 


or+0.08 
u - oo -0.09 


3 OO+O' 60 
°- uu -0.58 


10527 


138.7 


203.6 


195 


1.04 


250 


375 


8.0+°,;* 


30+ 009 
u - JU -o.io 


4-75i; 7 7 I 


7838 


262.6 


177.0 


161 


1.08 


375 


500 


12.5+1! 


4O+0.27 
u, ^ a -0.24 


i-86iS:!S 


5466 


378.0 


135.6 


121 


1.09 


500 


750 


16.9+Ji 


20+0' 45 
u - zu -0.20 


07+ 2 - 11 
u ' u '-0.07 


7003 


1123.4 


159.2 


169 


0.94 


750 


1250 


26 -°-ll 8 2 


o oo+o- 52 

u,uu -0.00 


00+ 029 
u ' uu -0.00 


8322 


3075.7 


198.0 


199 


1.00 


150 


250 


o 7+0.6 
a -'-0.6 


o sti+o- 09 


2.30 


10527 


138.7 


205.0 


196 


1.05 


250 


375 


io.o±° : » 


32+0 .12 
u "-> z _0.12 


2.30 


7838 


262.6 


188.1 


165 


1.14 


375 


500 


ll-9±i"S 


n «o+0.26 
U - 4S -0.23 


2.30 


5466 


378.0 


135.8 


125 


1.09 


500 


750 


13 1+ 21 


22+ - 21 
U - ZZ -0.12 


2.30 


7003 


1123.4 


162.3 


170 


0.95 


750 


1250 


i2-5±i;° 


n 9 o+0.36 
U - ZO -0.23 


2.30 


8322 


3075.7 


207.6 


200 


1.04 



Table 2. Model parameters for isothermal models using combined sets of data. The spectra from annuli around the cluster centre as 
used in Table hi have been fitted simultaneously using the same absorbing column density and metallicity for all annuli in the range 
ri — » r2 and by allowing for different temperatures in each of the annuli of Table hi The last three columns, as in Table hi describe the 
number of degrees of freedom and the quality of the fits. 



r\ (kpc) r2 (kpc) Z (solar) 



Nu 



DOF 
(10 20 cm- 2 ) (10 2 






50 


29+0- 05 
u - za -0.04 


50 


150 


31+ 006 
U - CSi -0.05 


150 


1250 


31+"- 07 
u -° -0.06 



3 94+0.51 
° ,y -0.51 
9 SS + 0- 44 
2 - SS -0.43 
-+0.38 



2.55^ 



0.36 



367.0 
422.0 
898.9 



285 
367 

857 



1.29 
1.15 
1.05 



mately constant with radius. We find tentative evidence that 
the isothermal model requires excess absorption above the 
Galactic value in the inner 100 kpc. 

As an aside, we mention that if we apply model (llj) 
to the spectrum extracted from a circle with a radius of 
30arcsec (150 kpc) and leave the Ni abundance as an ad- 
ditional free parameter, we find a Ni abundance of 2.3 ± 
0.8 times solar in this region. This relatively high Ni abun- 
dance is relevant for models of Type la supernov ae and 
their e nrichme nt of the intergalactic medium (e.g 
& Whhle 2000). 



Dupkc 



4.4 Spectral deprojection analysis 

The results discussed in the previous section are based on 
the analysis of projected spectra. In order to determine the 
effects of projection, we have also carried out a deprojection 
an alysis of the Chandra spectra using the method described 
by Allen et al. (2001a). This method decomposes the ob- 
served annular spectra (Table nj) into the contributions from 
the X-ray gas emission from nine spherical shells. 

The data for all nine annular spectra were fitted simul- 
taneously in order to correctly determine the parameter val- 
ues and confidence limits. The spectral model used therefore 



has 2 n+ 1 free parameters (where n — 9 is the number of an- 
nuli), corresponding to the temperature and emission mea- 
sure in each spherical shell and the overall emission- weighted 
metallicity (the metallicity is linked to take the same value 
at all radii, yielding a metallicity Z = 0.32+ 004 ). The ab- 
sorbing column density was fixed at the Galactic value. The 
temperature profile determined with the spectral deprojec- 
tion code is shown in Fig. pi 



5 X-RAY MASS ANALYSIS 

5.1 The mass model 

The observed X-ray surface brightness profile (Fig. H) and 
deprojected X-ray gas temperature profile (Fig. fel) may to- 
gether be used to determine the X-ray gas mass and total 
mass profiles in the cluster. For this analysis, we have used 
an enhanced version of the depro jection code described by 
White, Jones and Forman (1997[) and have followed a sim- 



ilar me thod to that described by Allen, Ettori fc Fabian 



(2001a ). Spherical symmetry and hydrostatic equilibrium 



are assumed. 

A variety of simple parameterizations for the cluster 
mass distribution were examined, to establish which could 
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Figure 7. As in Fig. 0, but for metallicity Z (in units of solar 
abundance) and photoelectric absorption column density JVh as 
determined from the combined fits in Table ti to X-ray spectra 
extracted from annuli at several radii. In the lower panel, the 
Galactic absorption value from H I studies is indicated by the 
dashed line. 
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Figure 8. Deprojected temperature profile of Abell 1835. The 
solid lin e is the best-fit temperature profile determined in 
Sect. |t[ 



provide an adequate description of the Chandra data. We 
find that the Chandra data can be modelled using a Navarro 



et al. (: 997) (NFW) mass model 



p(r) = 



Pcrit C 



(r/r s ) (1 +r/r s y 



(2) 



where p(r) is the mass density, p CI n = 3H(z) 2 /SirG is the 
critical density for closure of the universe (H(z) is the Hub- 
ble constant and G is the gravitational constant) and 



S c = 



200 



3 ln(l + c) - c/(l + c) ! 



(3) 



where c is the concentration parameter. The normalization 
of the mass profile may also be expressed in terms of an 
effective velocity dispersion, a = V50H(z)r s c. 

We step the mass model parameters r s and a through 
a regular grid of values. (The scale radius r s was stepped 
between 0.1 Mpc and 2.5 Mpc and the effective velocity dis- 
persion a between 500 km s^ 1 and 2500 km s" 1 on a grid of 
81x81 models.) Given the observed surface brightness pro- 
file and a particular parameterized mass model, the depro- 
jection code is used to predict a temperature profile of the 
X-ray gas (we use the median model temperature profile de- 
termined from 100 Monte-Carlo simulations) which is then 
compared for each mass model with the observed, depro- 
jected temperature profile (Fig. |8|). 

We have determined the goodness of fit of each simu- 
lated temperature profile by calculating the sum 



X 



E 



all radial bins 



J- obs ■*■ n 



(4) 



where T b s is the observed deprojected temperature profile 
shown in Fig. H and T mo dei is the model temperature profile, 
rebinned to the same spatial scale using an appropriate flux 
weighting. 

Fig. H shows the X-ray gas temperature profile im- 
plied by the best-fitting NFW mass model (x 2 = 11-9 for 
r s = 0.64t°;^ Mpc, a = 1275±^|g kms" 1 and c = 4.0±° ;^) 
overlaid on the deprojected temperature profile. The y 2 
value is mostly due to the differences between model and 
data close to the cluster centre, and may at least in part be 
attributed to the asymmetric substructure seen in Fig. ||. 

Fig. shows a contour plot of the x 2 values calculated 
for the different NFW models in the parameter space. The 
minimum x 2 value is marked by a plus sign. Three contours 
have been drawn at differences of A\ 2 — 2.3, 6.17 and 11.8, 
corresponding to formal confidence intervals of 68.3% (1 a), 
95.4% (2 a) and 99.73% (3a) for the two parameters (e.g., 



Press et al. 1992, page 661). Fig. llOl shows the 1 a range of 
enclosed masses for models around the best-fit model both 
for the total mass inside a sphere of radius r and for the total 
pro jected mass insid e a radius r (using the mass formula 
by Bartclmann 1996). The mass inside a sphere with the 
virial radius r2oo = c x r s for the best-fit model is M200 = 
2.7t° - 5 7 x lO^Meyr" 1 . 

We note that it is also possible to carry out this analysis 
with a nonsingular isothermal mass density profile 



p(r) = 



1 



2ttG r 2 +r 2 



(•>) 



where ai so is the velocity dispersion and r c is the core ra- 
dius. This yields parameters r c = 65.0lj ° kpc and o"i so = 
1280±60kms~\ The fit is slightly better, x 2 = 6.8, but the 
conclusions we can draw from both profiles regarding the 
total mass and X-ray gas profiles are very similar. 
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Figure 9. Contour plot of the 68.3% (Iff), 95.4% (2ff) and 99.73% 
(3<t) formal confidence contours on the NFW scale radius r s and 
effective velocity dispersion a. For this, a grid of 81 X 81 mass 
models was graded using the x 2 statistic by comparing the pre- 
dicted temperature profile with the deprojected Chandra temper- 
ature profile. Also shown for reference are the isocontours of the 
NFW concentration parameter c. The minimum x 2 model with 
X 2 = 11.9 is marked by a plus sign. 
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Figure 10. 1 a range and average value (dashed line) of the en- 
closed mass inside a sphere of radius r and the total projected 
mass inside a radius r. 



5.2 Lensing predictions from the X-ray mass 
models 

In Fig. hll (left panel) an optical B-R colour image is shown 
of the central 1.5x1.5 arcmin of Abell 1835. The image has 
been constructed from the 20min B-band image shown in 
Fig. ti and a lOmin R-band image from the same observing 
run that was also taken from the CFHT archive. In this 
greyscale image, blue objects appear black, and red objects 
appear white. The contrast has been stretched that very 
faint objects can be seen. The two vertical stripes are due 
to bleeding columns in the R-band image that fortunately 
did not affect the cD galaxy. 

In the central part of Fig. [Ill (left panel) , it can be seen 
that only the cD galaxy appears blue. The two arcs to the 
south east and south west of the cD galaxy were already 
pointed out in Fig. 0. The larger arc to the south east was 
discovered by Edge et al. (in preparation). The other arc 
has not been discussed in the literature so far. No redshift 
has been determined for either object in the literature, but 
their optical appearance and colours indicate that they are 
gravitationally lensed images of background galaxies that 
are situated at some distance be hind the galaxy cluster (e.g., 
Schneider, Ehlers fc Falco 1992; ). 

In the R-band image (not shown) we find a further pair 
of very thin and elongated features parallel to, but 9 arcsec 
south of the south east arc (termed CI and C2 in the follow- 
ing) . This pair of arcs is too faint to be seen in Fig. hll (left 
panel). The total length of the pair is comparable to arc A, 
but is only barely visible in the B-band image. If we assume 
that the very faint B-band brightness of the pair is due to 
the redshifting of the Lyman break out of the B-band filter 
used (3800-4800 A), w e may deduce a red shift of z > 4.0 for 



this object pair (e.g., Madau et al. 1996) 



The separations of the centres of arcs A and B from 
the centre of the cD galaxy are 8a ~ 31.1 arcsec and Ob ~ 
25.0 arcsec. In a circularly symmetric gravitational lens, arcs 
appear close to the Einstein radius #e 



-Dclu 



4GM 



(6) 



Schneider et al. 1992) 



-^clus -L'arc l- 

where M is the enclosed mass (e.g.. 
D C \ UB = 1.02 Gpc, D arc and -D c ius-arc are the angular size 
distances of the cluster, the arc, and between cluster and arc. 
G is the gravitational consta nt and c is the spe ed of light. 
For the case of Abell 1835, |Allen et al. (1996[) calculated 



the mass in a circular aperture inside arc A (their Fig. 13) 
and found a value between 1.4x 10 14 M and 2x 10 14 M for 
values of the redshift of the arc 0.6 < z a rc < 3.0. Similar 
values follow for arc B. 

The NFW mass models in the region of best goodness 
of fit in Fig. have only little variatio n of the total mas s 
inside 31.1 arcsec (153 kpc): If we use the Bartelmann (199C) 



formulae for the total mass in a circular NFW profile, the 
lcr formal confidence interval data ranges from Mx = 1-06 x 
1O 14 M to M x = 1.27 x 1O 14 M (Fig. [HJ. The values 
agree wi th the 1.26 + n ? 2 x 10 14 M determined using ROSAT 



data by Allen (1998). The values are lower, however, than 
the estimate from the lensing argument assuming circular 
symmetry. 

From the X-ray image in Fig. H we can deduce that in 
the inner parts the X-ray isocontours are elliptical with an 
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Figure 11. (left panel) B-R colour difference image of the central 1.5xl.5arcmin of Abcll 1835 from archival CFHT data. North is 
up and East is to the left. The greyscale-contrast has been stretched to show blue objects in black and red objects in white. The two 
vertical features to the left and right of the central cD galaxy are due to bleeding CCD columns in the R-band image that were caused 
by bright stars to the north of the field. We have replaced the pixel values in these columns by the median background value. Features 
of particular interest in this image are the very blue central part of the cD galaxy and the two arcs to the south east (A) and south west 
(B) of the galaxy. We also note that to the north of the cD galaxy a straight blue feature is seen (arrow), (right panel) Contour plot of 
gravitationally lensed images by an elliptical NFW profile with a concentration parameter c = 3.3, scale radius r s = 900 kpc, ellipticity 
e = 0.2 and position angle of 330 degrees (measured counterclockwise from north). The cluster centre is marked by the cross. See the text 
for details on source positions, redshifts and sizes. The images have been labelled with upper case letters. B' is the predicted counter arc 
of arc B. 



axis ratio of / = 0.85 (Sect. pt). This suggests that also the 
mass distribution of Abell 1835 is not circular ly symmetric. 
Numerical simulations by Bartelmann (1995) showed that 
the lensing mass estimate under the assumption of circular 
symmetry needs to be corrected towards lower values. As an 
average correction factor he quotes 1.6, but this may vary 
downwards to unity and upwards to ~2 from case to case. 



with constants as described after eq. (pi). One can then gen- 
eralize eq. (M) for a mass distribution with elliptical symme- 
try by making the transformation 



+ 



For the lensing mass estimate the ellipticity thus could be 
an important factor, in order to go beyond the approxima- 
tion of circular symmetry, we need to use an elliptical mass 
model. 

The normalized surface mass density k = S/ S cr j t of 
a spher ical NFW mass distribution was derived by Bartel- 
mann (1996). Using the ratio x — r/r s of radius r and scale 



(l + c)a (1 
(Schramm 1990J, |l994t JKassiola fc Kovner 1993 



(10) 



Schneider and Bartelmann 1994) where xi and x-2 are carte- 



Kormann 



radius r s this is given by 



k(x) = 2k s 



where 



/(*) = 



fix) 







v^ 



2= arctan , /^-i ( x > 1) 

.2_j y se+1 v ' 

arctanh i/yxf ( x < 1) 



(7) 



(8) 



(x = l) 



(k (1) = I « s ) and k b = S c pcrit r B £ cr \ t ■ The critical surface 
mass density E cr it is given by 



2-icrit — 



Da 



4 TV G Ddus Adus 



(9) 



sian coordinates and e is the elliptical parameter. The ellip- 
tical parameter e is related to the axis rati o / via / = y^ . 
Using the summation method developed by Schramm (1994) 
we can then construct gravitational lens models for the ellip- 
tical NFW mass model within the viable range of parameters 
from Fig. H. We use an axis ratio / = 0.8 and a position an- 
gle of 330 degrees (measured counterclockwise from north) 
taken from the optical appearance of the cD galaxy. We do 
not use the axis ratio of the X-ray gas because the gas pres- 
sure causes the gas distribu tion to be more circular than 
the mass distribution (e.g., puote fc Canizares 1996 ). We 
fix the centre of the mass model to the optical centre of the 
cD galaxy. 

Without redshift information about the arc candidates, 
any parameterized lens modelling of the system is degener- 
ate with respect to the mass normalization of the cluster. 
It is also necessary to make an assumption about whether 
to interpret the arcs as single images which have been elon- 
gated by the lens effect, or as fold-arc s consisting of two mir - 
ror images of the same object (e.g., Schneider et al. 1992 ). 
Since we do not detect counter images for arc A, it is likely 
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that it is a single, elongated image of background galaxy. 
We also cannot detect a counter image for the pair C1/C2 
in the CFHT R-frame, which could indicate that they are 
bright spots in a common, darker background galaxy that 
appears elongated by the lens effect. On the basis of the 
archival CFHT data, however, it is impossible to identify 
unambiguously which images belong together, or to rule out 
the existence of further counter images. 

We have first studied the lens effect of an NFW model 
with a = 1460 km s" 1 and r B = 0.9 Mpc, which is situated 
inside the 1 a confidence region in Fig. bl With this mass 
distribution, arc A can be explained as a single elongated 
image of a background galaxy at z\ ~ 2.7. The arc pair 
C1/C2 in this model are images of two objects at zq ~ 5.0 
with a projected physical (unlensed) separation of ~ 10 kpc. 
This redshift is consistent with our limit zc =S 4.0 which 
was obtained from the fact that the pair is not visible in the 
CFHT B-frame. An arc at the position of arc B is predicted 
for several source redshifts, but the most interesting inter- 
pretation is a background object at zb ~ 2.25, because at 
this redshift the lens model predicts counter images in the 
form of a radial arc at the position of the straight blue fea- 
ture to the north of the cD (Fig. |llj, left panel; t he existence 
of r adial arcs in NFW profiles was shown by Bartclmann 
199dj~" 



The image configuration predicted by this model is 
shown in the right panel of Fig. O. Circular Gaussian bright- 
ness profiles have been assumed for all sources with a full 
width at half maximum of 0.45 arcsec for source A and 
0.25 arcsec for sources B, CI and C2. The contours have 
been spaced at 10%, 30% and 50% of the central surface 
brightness. The assumed source positions on the sky with re- 
spect to the cluster centre in the form (Ax, Ay) are: source 
A: (-5.06", -6.20"), source B: (0.50", -2.44"), source CI: 
(-7.05", -9.03"), source C2: (-5.78", -10.22"). 

Secondly, since arc A appears roughly symmetric 
(Fig. O, left panel), we have also investigated the possibility 
that it is a fold-arc. This interpretation may be supported 
by the similarity of arcs CI and C2. We find, however, that 
none of the models inside the 1 a confidence region in Fig. H 
are massive enough to explain a fold-arc at the radius of arc 
A (31.1 arcsec or 153 kpc). In this case it would be necessary 
to scale the surface mass density n —* k — q x k of all mod- 
els in the formal 2 a contour by correction factors between 
q = 1.1 and q = 1.4. The redshift of arc A in such models 
would typically be z\ ~ 2.0. (The redshift is constrained to 
be around za ~ 2.0 because it is also necessary to explain 
arcs B and C1/C2. Arc B would have a redshift zb ~ 0.9. 
CI and C2 would be interpreted as images of a source at 
zc > 4.2. The radial feature to the north of the cD could 
again be interpreted in this model as a radial arc.) 

We have also studied different axis ratios within the 
context of the fold-arc model: Lowering the axis ratio re- 
duces the correction factor q because the ellipticity causes 
the critical curve to extend further out. We thus find that 
elliptical / ~ 0.6 models within the 2<r confidence contour 
of Fig. are able to explain the CFHT data. The axis ratio 
should be larger than / = 0.5, however, because otherwise 
it is not possible to have an elongated arc at the position 
of arc B (without additional perturbations of the potential). 
By using nonsingular isothermal mass profiles (eq. (ph) the 
correction factor q can also be reduced further. 



Our analysis clearly favours the first, single-image, in- 
terpretation of arc A, as well as C1/C2, because the fold- 
arc model predicts counter-images for arcs A, C1/C2 that 
are not observed. The available strong lensing data are ex- 
plained by a mass model entirely consistent with the Chan- 
dra temperature profile, indicating that there are no sig- 
nificant contributions from non-thermal sources of pressure 
(e.g., bulk motions, magnetic fields). The fold-arc models 
require special tuning of the mass distribution, although it 
may be possible to explain the additional mass with sys- 
tematic uncertainties of our X-ray mass analysis due to the 
residual substructure seen in the X-ray emission, or by a 
merging subcluster or some other line-of-sight enhancement 
present in this region. Nevertheless, the predicted counter 
images would have to be found before this option can be 
considered further. 

We have not attempted any more detailed lens mod- 
elling including sub-halos due to cluster member galaxies 
because the available data are of insufficient quality to un- 
ambiguously establish the arc properties and identify image 
pairs. A more rigorous analysis aiming to place precise con- 
straints on the redshifts of the arcs and corresponding arc 
properties is a subject for future st udy and will best be done 



with the Hubble Space Telescope (Kneib 2000; smith et al 



200C). Moreover, spectroscopy will be required in order to 
ensure that th e arc candidates are indeed galaxies at high 
redshift (e.g., Ebbels et al. 1998). A measurement of the 



redshift of arc B alone would provide a very strong test of 
the mass models discussed here. 



5.3 The X-ray gas mass fraction 

The X-ray gas-to-total-mass ratio as a function of radius, 
/gasM, determined from the Chandra data is shown in 
Fig. [12. For this figure, the gas mass fraction was calcu- 
lated for all models contained in the 1 a confidence contour 
of Fig. H. The dashed line corresponds to the best-fit model, 
the solid lines are the envelopes of all calculated profiles. 
The / gas value rises rapidly with increasing radius within 
the central ~ 50 kpc and then flattens, remaining approx- 
imately constant out to the limits of the data at r — 1.0 



Mpc, where we measure / g: 



0.171] 



U.U^U r — 1. 
0.022 a 50 



Following the usual arguments, which assume that the 
properties of clusters pro vide a fair sample of those of the 
Universe as a whole (e.g., [White, Efstathiou fc Frenk 19931; 



White & Fabian 199E 



Evrard 1997 



Ettori fe Fabian 1999) 



we may use our result on the X-ray gas mass fraction in Abell 
1835 to estimate the total matter density in the Universe, 
fl m . Assuming that the luminous baryonic mass in galaxies 
in Abell 18 35 is approximate l y one fifth of the X-ray gas 
mass (e.g., White et al. 1993, Fukugita, Hogan & Peebles 



199£ ) and neglecting other possible sources of baryonic dark 
matter in the cluster, we obtain Q m = (fib/l-2/gas), where 
fib is mean baryon density in the Universe and /150 is the 
Hubble constant in units of 50kms -1 Mpc -1 . For fib^io = 
0.0820 ± 0.0072 ( JQ'Meara et al. 2000| ), we measure fi m = 
0.40±0.08/i7 n - 5 . 



THE PROPERTIES OF THE COOLING 
FLOW 
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Figure 12. Upper and lower 1 cr limits on the gas mass fraction 
fraction as a function of radius r. In addition, for reference the 
gas mass fraction has been plotted also for the best-fit model 
r s = 0.64 Mpc, a = 1275 km s -1 (dashed line). 




radius 



Figure 13. Electron density n e and cooling time t coo i as a func- 
tion of radius r. The cooling radius r coo i at which the cooling time 
profile becomes larger than the age of the universe | Ho = 13 Gyr 
(dashed line) is r coo i ~ 230 kpc. 



6.1 Density and Cooling times 

Using the density and temperature pro files determined with 
the best- fit NFW mass model in Sect. 



5.1 



we can calculate 
the cooling time of the intracluster gas (i.e. the time taken 
for the gas to cool completely) as a function of radius. The 
results on the electron density and cooling time profiles are 
shown in Fig. [13. The overall profile of the density data can 
be approximated by a /3-model 



n c (r) oc 



1 + 



(0 



-, _i£ 



(11) 



6.2 Standard cooling flow modelling and image 
analysis 

In this section we report on the results from fitting cooling 
flow models to the spectra extracted from full circles centred 
on the geometrical cluster centre in the 0.5-7.0 keV band. In 
the standard cooling flow picture this yields the integrated 
mass deposition rate of the matter that cools out of the flow 
within a given radius. 

We model the cooling flow using the RJCOOL e mission 
model constructe d according to the prescription in lohn- 



with a bore radius, rv = 39 ± 3 kpc. B = 0.53 ± 0.01 and a stone et al. (1992J ) and include intrinsic absorption that may 



central density, n t 

Trip Hashed 



(0) = 10.2 ±0.8 x 10~^ cm^. 

inp in the lower panel nf Fig [l 3J indr- 



be acting on the cooling flow in the core of the cluster with a 



cates t he age of the universe for our chosen cosmology with 

Ho = 50 km a -1 Mpc -1 For the assumed Galactic column 



ZPHABS absorbing model (Balucinska-Church & McCam- 



densitv l of 2.3 x 10 20 cm" 2 ([Dickey fc Lockman 19!)(|) 



mon 1992| ) at the redshift of the cluster. We use the MEKAL 
plasma emission model (Kaastra & Mewe 1993; 



Liedahl et 



al. 1995) to model the gravitational work done on the gas as 



measure a central cooling time of t C ool ~ 3 x10 s yr, and a 
cooling radius, at which the cooling time first exceeds the 
age of the universe, of r coo i ~ 230 kpc. In Fig. hj we plot the 
entropy S = — |W as a function of radius, calculated using 

the temperature and electron density profiles of the best-fit 
model. 

The fact that cooling times in the centres of cooling flow 

clnsters, are comparable to nr even significantly less than the 



well as emission of the intracluster gas outside the cooling 
flow: 



MODEL 2 



PHABS ga i x 



MEKAL (T U gh; Z; K ) 



+ ZPHABS (AJV H ) x RJCOOL (M)] .(12) 



'S_£ 

age of t | he universe is well known _i^ea^^al^l^7^L Silk 1976 
Cowie & Binney 1977; Fabian & Nulsen 1977) and was an 



The free parameters are given in brackets. The ove rall ab- 
sorption in this model is fixed at the Galactic value (Dickey 



Lockman 1990). The intrinsic absorption ATVh acts only 



important motivation for the development of cooling flow 
models. In the next sections we explore the properties of 
such a cooling flow model for Abell 1835. 



on the cooling flow emission. The normalization of the cool- 
ing flow model, M, measures the rate of gas cooling out of 
the flow. The other parameters of the cooling flow model, 
namely the metallicity Z and upper temperature Thigh, are 
fixed to the values of the MEKAL com ponent. This model 
is identical to model C in Allen (200C). It is a multi-phase 
model in the sense that at each radius all temperatures from 
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Figure 14. Entropy S as a function of radius. To guide the 
eye the entropy value for a temperature of lOkeV according to 
the numerical relation from Ponman, Cannon & Navarro (1999) 
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has been plotted with dash ed lines: S=567.0 at 0.1 virial radii, 
0.1 cr s = 257kpc (Sect. O). 



Figure 15. Integrated mass deposition rate M as a function of 
radius r determined spectrally (open circles) and from the image 
analysis (filled circles). The solid line is a broken power-law fit 
to the image analysis profile with a break at rij rca k = 42 kpc and 
exponents of 2.4 before and 0.4 after the break. 



the upper temperature down to the lowest temperature are 
present. 

In Table [3| the results of the model fits are given for 
outer radii from 5arcsec (25 kpc) to 50arcsec (250 kpc). In 
Fig. |l_5J we have plotted with open circles the mass depo- 
sition results with an even finer spacing of the extraction 
circles. All of the fits are statistically acceptable, with re- 
duced x 2 values around or slightly above unity. There is no 
significant improvement in \ over the simpler isothermal 
models eq. (hi). The Chandra data for Abell 1835 cannot 
discriminate between single phase and multi-phase models. 

Also shown in Fig. hq (solid circles) is the mass depo- 
sition rate that may be deduced from t he X-ray luminos- 
ity distribution (e.g., White et al. 1997 ) using our best-fit 
mass model from Sect. 



3.1 



It can be seen that the results 
from the spectral modelling and from the image analysis 
agree rather well out to 30 kpc, where the spectral mass 
deposition profile has a break at a mass deposition rate of 
231I53 Mq yr _1 . The profile from the image analysis then 
itself has a break around 40 kpc. The solid line corresponds 
to a broken power-law fit to the mass deposition profile from 
the image analysis. Both profiles rise with a similar slope af- 
ter their respective breaks, but disagree with respect to the 
amplitude. The spectral mass deposition profile rises up to 
a value of ~ 500 Mq yr _1 within 250 kpc, whereas the pro- 
file from the image analysis rises to values of 12 50 Mp, vr~ ' 



which is consistent with the value from ROSAT ( Allen 2000 ) . 
The intrinsic absorption AJVh tabulated in Table H act- 
ing on the cooling flow emission rises from a mean equivalent 
hydrogen column density of AA<h ~ 0.2 x 10 22 cm -2 at a 
radius of 50arcsec (250 kpc) by a factor of ~ 1.5 when we 
look deeper into the cluster core. The statistical error bars 
on the intrinsic absorption are large, however. The tempera- 



ture of the MEKAL component (which is linked to the upper 
temperature of the cooling gas) rises from about 5 keV in the 
centre to 9 keV within a circular aperture of 750 kpc radius. 

6.3 Interpretation 

6.3.1 A young central cooling flow? 

Standard co moving multiphase cooling flow models (e.g., 
Nulsen 1986) typically predict the profile of the mass de- 
position rate as a function of radius to be a power law. A 
break in the profile of the mass deposition rate, as seen in 
Fig. hsl may suggest tha t the flow mechanis m is changing 
around the break radius. Allen et al. (2001b) interpret the 



presence of such breaks as an indicator of the age of the cool- 
ing flow, which they suggest is likely to be closely related to 
the time since the central regions of the cluster settled down 
following the last significant disruption by merger events or 
AGN activity. Given the cooling time at a radius of 30 kpc 
(Fig. [13J) , this would then limit the age of the cooling flow 
in Abell 1835 to be < 6x 10 8 years. 

The natural question which must be asked is whether 
such a young age for the cooling flow in Abell 1835 is plausi- 
ble. Given the clear substructure observed in the core of the 
cluster (Sect. H), some disturbance in the recent past seems 
likely. As discussed in Sect. H, the small emission front seen 
in Fig. to the south-east of the cluster core may have 
be en caused by a cool, inf alling gas cloud (e.g., see the work 
of iMarkevitch et al. 2000|). A 31.3 mjy FIRST radio source 



(Becker et al. 200C) is located in the cD galaxy. Although 



this source is relatively quiescent at present, it may have 
been much more powerful in the recent past and injected 
energy and momentum into the intracluster medium (pos- 
sible problems with radio source heating are discussed by 



12 Schmidt, Allen & Fabian 



Table 3. Cooling flow model fits (eq. (112)) to spectra from circular regions with radii r. We list the best-fit values for the model 
parameters mass deposition rate M, upper temperature Thigh, metallicity Z and intrinsic equivalent hydrogen column density ATVjj. We 
also list the x 2 value, the number of degrees of freedom (DOF) and the reduced \ 2 value (see Table 
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David et al. 2001 and Fabian et al. 2001). One should also 



recall that within the context of structure formation mod- 
els, a cluster as hot and massive as Abell 1835, observed at 
a redshift z ~ 0.25, is highly unlikely to have formed more 
than a few Gyr before. 

If the cluster remains undisturbed, the results on the 
cooling times and the mass deposition rate in Figs. Il3| and ha 
suggest that the core of the cluster will develop into a very 
massive cooling flow within a few Gyr. It is also possible 
that the cluster contained a large cooling flow before its 
most recent disturbance. Note that similar studies of the 
cooling flows in the Abell 2390 and 3C295 clusters indicate 
older ages of 2-3 and 1-2 Gyr, respectively. 

We note that the fact that the break appears at a 
slightly smaller radius in the spectral analysis than the im- 
age analysis could simply be a consequence of the fact that 
although the gas starts to cool appreciably around 40 kpc, 
it has only had sufficient time to cool completely within 
r ~ 30 kpc. Note also that given the complex central mor- 
phology in Fig. 0, simple, spherically-symmetric cooling-flow 
models are unlikely to be applicable in detail. 

As discussed by Allen et al. (2001c), in cases where the 



temperature of the cluster gas continues to rise beyond the 
outer edge of the cooling flow, spectral-measurements of the 
mass deposition rate covering regions larger than the cool- 
ing flow and using models similar to eq. (|12[), are likely to 



overest: mate the true mass deposition rate. This is because 



the presence of relatively cool, ambient spectral components 
will tend to be modelled as part of a cooling flow. This sit- 
uation is likely to have affected previous ASC A studies of 



the integrated cl uster spectrum for Abell 1835 (Allen et al 



1996 ; Allen 2000). We also note that the application of such 
models to large spatial regions may lead to the detection of 
low-temperature 'cut-offs' in the cooling flow spectra, which 
would simply be due to the fact that complete cooling occurs 
only within the innermost cooling flow proper. 



6.3.2 XMM observations of Abell 1835 



In a recent paper Peterson et al. (2001) present new re- 
sults on grating spectra of Abell 1835 taken with the Reflec- 
tion Grating Spectrometer (RGS) aboard the XMM-Newton 
X-ray observatory. They find no evidence for line emission 
from gas below 2.7 keV and, neglecting intrinsic absorption, 
place a 90 per cent confidence upper limit on the mass 
deposition rate from the cooling flow of 200 M0yr -1 for 
Ho = 70kms _1 Mpc -1 , On = 0.3 and A = 0.7. In fact, 
RGS has not found convincing evidence for a cooling flow 
in any galaxy cluster yet. Due to the inverse square depen- 



dence of the mass deposition rate on the luminosity dis- 
tance, the RGS limit translates into 315 MQyr -1 for the 
Ho — SOkms -1 Mpc -1 Einstein-de Sitter cosmology used 
here. 

Overall, the RGS result for Abell 1835 implies a range 
of temperatures from 9 to 3 keV within the inner 0.5 ar- 
cmin radius (150 kpc), with the emission measures of that 
gas roughly matching a cooling flow. This is in agreement 
with our findings which show a radial gradient in the am- 
bient gas temperature from ~8keV at 150 kpc to ~4keV in 
the cluster core. In addition, we find evidence for complete 
cooling within the inner 30 kpc. 

The RGS limit on the mass deposition is not easily rec- 
onciled with the large (M > 1000 Mq yr _1 ) mass deposi- 
tion rates inferred from previous studies of the integrated 
cluster spectrum observed with ASCA and ROSAT. Within 
the context of the young cooling flow model discussed above, 
however, the RGS results can be more easily accommodated. 
If one also accounts for the effects of intrinsic absorption, 
as required by the Chandra data, the apparent RGS limit 
of 315 Mq yr _1 can be increased and is consistent with the 
Chandra result of 231+^ M Q yr" 1 . 

We note that the revised mass deposition rate of ~ 
200-300 M yr" 1 within r ~ 30 kpc in Abell 1835 is compa- 
rable to the optical star formation rate (normalized to the 



infer r ed number of A star s within the same region; Allen 
1995; Crawford et al. 1999]). Thus, it appears that a signifi- 



cant fraction of the material being deposited by the cooling 
flow within the central 30 kpc of Abell 1835 forms stars with 
a relatively normal IMF. pdge (2001 ) also reports the de- 
tection of a substantial molecular gas mass (~ 2 x 10 x Mq) 
within the central regions of the cluster. 



6.3.3 The origin of the temperature profile 

The temperature profile in a cluster core is a result of the 
thermodynamic history of the intracluster gas. Chandra en- 
ables us for the first time to compare in detail the ob- 
served temperature profiles in clusters with those pred icted 
by numerical studies o f galaxy cluster fo r mation (e.g., |Eke, 
Navarro fc Frenk 19 98J |Frenk et al. 1999|; [Lewis et al. 200' " 
2000D 



Pearcc et al. 2000[). I n the studies by |Eke et al. (199., , 
and Frenk et al. (1999), the temperature profiles of massive 
galaxy clusters (determined with smooth particle hydrody- 
namics, without including the effects of cooling), are often 
observed to drop slightly towards to core (r < 100 kpc), 
then remain roughly constant on a plateau to r ~ 1 Mpc, 
before dropping again at larger radii. 
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Pearce et al. (2000) present temperature profiles for a 
sample of 20 simulated clusters, both with and without the 
effects of radiative cooling. For their largest cluster, the re- 
sults obtained wit hout including t he e ffects of cooling ar e 



similar to those of [Eke et al. (1998J ) and |Frenket al. (1999| ) 
However, when cooling is included, these authors find a 
"precipitous decline in temperature" in their largest cluster 
which they deem evidence for a cooling flow. The tempera- 
ture profile for this cluster appears similar to that observed 
for Abell 1835 (Fig. |6|). The work of Pearce et al. therefore 
suggests that the sharp temp erature drops in the cores of 
cooling flow clusters (see also Allen et al. 2001a; David et 
al. 20011) could be due to the effects of radiative cooling over 
the lifetime of the cluster. Since in other simulations with 
cooling (e.g., Lewis et al. 200C) the temperature profile 



seems to keep rising with decreasing radius we note, how- 
ever, that the effect of cooling in numerical simulations on 
the cluster temperature profiles does not seem to be entirely 
clear yet. 



6.3.4 A massive cooling flow? 

In this section we address whether an old (~5 Gyr) , very 
massive (M ~ 1000 Moyr -1 ) cooling flow, encompassing 
the entire central 150 kpc radius, can e xist in Abell 1835. 
In this case, as mentioned in Sect. 3.3.2, our problem is to 



reconcile the large mass flow rate with the relativ ely small 



amount of X -ray l ine emission observe d with RGS. Peterson 
et al. (2001 ) and Fabian et al. (2001 ) discuss a number of 
processes that may be relevant, two of which are investigated 
here. 

The first model invokes an inhomogeneous distribution 
of metals that can be entered into XSPEC using the follow- 
ing prescription: 

MODEL3 = PHABS g ai x (13) 

[MEKAL (T h i gh ; Z; K) + ZPHABS (AJVh) x 
(RJCOOLi (M; Z = 0) + RJCOOL2 (0.1 M; 10xZ))] 

where the free parameters (given in brackets) are the same 
as in eq. (112), but with the mass deposition rate distributed 
between two cooling flow components of which one has zero 
metallicity and the other all the heavy elements. This model 
describes a situation where most of the heavy elements are 
confined to small regions, so that the bulk of the cooling flow 
is extremely metal poor. The ratio of the mass of the metal 
poor and metal rich regions was fixed as 10:1. We have fitted 
this model to the Chandra spectrum extracted from a circle 
with a radius of 30arcsec (150 kpc). The best-fit parame- 
ters are given in Table Id. The total mass deposition rate is 
1070 A/0 yr _1 . Although the mass deposition rate is large, 
however, the emission lines due to cool gas are suppressed. 
The FeXVII line at 15. 01 A, for example, is suppressed by 
about a factor 4. However, the Fe XXII and Fe XXIII com- 
plex around 12.2 A, which is strongly constrained by the 
RGS data, is only suppressed by a factor 2 by this model. 
This appears to be inconsistent with the RGS result. 

The second model incorporates absorption that only 
acts on the gas below a certain temperature threshold and 
can be entered into XSPEC using the following prescription: 



[MEKAL (Thi gh ; Z; K) + RJCOOLj (M; Ti ow -*T higll ) 4 

ZPHABS (AiV H ) x (RJCOOL2 (M;0keV->Ti ow ))] . 

Again, the model is very similar to the one presented in 
eq. ([12|). The best-fit values for the spectrum extracted from 
a circle with a radius of 30arcsec (150 kpc) are also given 
in Table W. The emission measure of gas below 2.7 keV is 
suppressed in the same way as the previous model, but this 
model also fails to adequately suppress the Fe XXII and Fe 
XXIII emission around 12.2 A. 

Although both models have significantly reduced emis- 
sion from emission lines by cool (< 2.7 keV) gas, their emis- 
sion measure from gas at intermediate temperatures does 
not seem to be consistent with the RGS results. 



7 CONCLUSIONS 

We have presented the analysis of 30 ksec of Chandra ACIS- 
S3 observations of the galaxy cluster Abell 1835, of which 
20 ksec were usable for spectroscopic analysis. Overall, the 
X-ray image shows a relaxed looking galaxy cluster. The X- 
ray isocontours in the centre are elliptical with an axis ratio 
of 0.85. In the inner 30 kpc radius, however, we detect clear 
substructure in the image. 

Using isothermal fits to spectra from annular regions 
around the cluster centre we have measured the tempera- 
ture, metallicity and photoelectric absorption of the cluster 
gas as a function of distance from the cluster centre. The 
temperature analysis reveals a steep drop in the tempera- 
ture from 12keV in the outer regions to only 4keV in the 
cluster core. We also find evidence for a rise in the photoelec- 
tric absorption towards the cluster core, but no indication 
of a change in the metallicity with radius. 

The Chandra data provide tight constraints on the grav- 
itational potential of the cluster which can be parameterized 
by an NFW mass model. The best-fit NFW model has a 
scale radius r B = 0.64+g'ig Mpc, an effective velocity dis- 
persion a — 1275 ± 160kms -1 , a concentration parame- 
ter c = 4.0±o'.64 and a total mass M = 9xl0 14 M Q within 
1 Mpc. We measure the X-ray gas mass fraction of Abell 
1835 as a function of radius and use this to deduce a cosmic 
matter density f2 m = 0.40 ± 0.09 /%, ' 5 . 

The projected mass within a radius of ~150 kpc implied 
by the presence of gravitationally lensed arcs in the cluster 
is in good agreement with the mass models preferred by the 
Chandra data. This rules out any significant contributions 
from non-thermal sources of X-ray gas pressure (e.g., bulk 
motions, magnetic fields). 

The detection of the Sunyaev-Zeldovich (SZ) effect 
( [gunyaev fe Zeldovich 1972) in A bell 1835 was recently re- 
ported by Mauskopf et al. (2000| ). The temperature (Fig. ||) 
and electron density (Fig. |13|, upper panel) profiles deter- 
mined in this paper can be used to predict the cosmology- 
dependent SZ effect in the cluster. This calculation and the 



cosmologica l implications of the SZ-detection by Mauskopf 
et al. (200C) are described elsewhere (Schmidt & Allen 2001, 



MODEL4 = PHABSgai x 



(14) 



in preparation). 

The Chandra data imply a radiative cooling time for 
the gas in the centre of Abell 1835 of about 3xl0 8 yr. By 
fitting cooling flow models to spectra from circular regions 
around the cluster centre, we measure a mass deposition pro- 
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Table 4. Best-fit parameters for possible massive cooling flow models. The upper limit on the metallicity in the first model was fixed to 
the software limit Z = 5 in the high metallicity component. 



model 



Thigh (keV) T low (keV) Z (solar) AN H (10 2 



2 ) M(M yr" 1 ) X 2 



DOF 



inhomogeneous mctallicities 
absorbed cool gas 9.9 



n.o+ti 

+7.1 
2.0 



2.6 



+0.7 
0.4 



0.50- 
0.29 4 



0.08 
0.51 



+ 0.02 

-0.01 

+ 0.1.') 

0.13 



975 +170 -° 
y<o.u_ 14g 

935.0±*&° 



322.3 
331.7 



285 
284 



1.13 
1.16 



file that can be represented by a broken power law with a 
break at a radius of 30 kpc, within which the integrated mass 
deposition rate is 23llgg Mq yr -1 . We observe a similar be- 
haviour in the mass deposition profile that can be derived 
from an image deprojection analysis, although the break 
in this case appears at a slightly larger radius r ~40 kpc. 
The spectral and the imaging mass deposition profiles agree 
within r ~30kpc. If we associate the cooling time of the X- 
r ay gas at the brea k radius with the age of the cooling flow 
( Allen et al. 2001b| ), we obtain an age of ~ 6x 10 8 yr. 

Recent results from th e XMM-Newton X-ray observa- 
tory (Peterson et al. 2001) limit the mass deposition rate 
in Abell 1835 to 315MQyr _1 , neglecting intrinsic absorp- 
tion. This limit appears to be consistent with our result of 
having a young cooling flow with a mass deposition rate 
M = 231^53 M©yr _1 in the centre of Abell 1835, particu- 
larly if one also includes the effects of intrinsic absorption as 
required by the Chandra data. The young cooling flow sce- 
nario is supported by the detection of substructure in the 
core, which shows that this region has not completely re- 
laxed. Strong signs of star formation have been detect ed in 
the inner 30 kpc with optical spectroscopy (Allen 1995) that 
could be fuelled by the inflowing gas. 
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